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Abstract. - We consider the problem of equilibration of a single oscillator system with dynamics 
given by the generalized classical Langevin equation. It is well-known that this dynamics can 
be obtained if one considers a model where the single oscillator is coupled to an infinite bath of 
harmonic oscillators which are initially in equilibrium. Using this equivalence we first determine 
the conditions necessary for equilibration for the case when the system potential is harmonic. We 
then give an example with a particular bath where we show that, even for parameter values where 
the harmonic case always equilibrates, with any finite amount of nonlinearity the system does not 
equilibrate for arbitrary initial conditions. We understand this as a consequence of the formation 
of nonlinear localized excitations similar to the discrete breather modes in nonlinear lattices. 
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[^ — . ' Introduction. — One of the simplest phenomenolog- 
T-H ical ways of modeling the interactions of a system with a 
'^heat bath is through the Langevin equation which, for a 



FD relation 



{Tl{t)fl{t')) = ksT^it - t') 



(3) 



\o 



single particle, of unit mass and moving in a one dimen- 
sional potential V{x), is given by 



C^ 



X — ~dV{x)/dx — jx + r]{t) , 



(1) 
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where 77 (t) is a Gaussian white noise which satisfies 
the fluctuation-dissipation (FD) relation {rj(t)rj{t')) — 
"IkBT^Slt — t') . Here T is the temperature of the heat 
bath. The dynamics in eq. (1) ensures that at long times 
the system reaches thermal equilibrium. Thus the phase 
space density P{x,p,t), where p — x, converges in the 
limit t — > 00 to the Boltzmann distribution e^"^'/Z 
where Hg = p^ /2 + V{x) and Z is the corresponding parti- 
tion function. The proof for this uses the correspondence 
between the Langevin equation and the Fokker- Planck 
equation [1]. 

However the J-correlated nature of the noise is un- 
physical and this has led to the study of the generalized 
Langevin equation [2,3] 



X = —dV{x)/dx 



dt'j{t - t')x{t') + ■q{t) 



(2) 



where the noise is correlated and the the dissipative term 
involves a memory kernel. The noise is again Gaussian and 
is related to the dissipative term through the generalized 



A standard method of microscopically modeling a heat 
bath is to consider an infinite collection of harmonic os- 
cillators, with a distribution of frequencies, coupled lin- 
early to the system. In that case it can be shown [4, 5] 
that the effective equation of motion of the system is pre- 
cisely given by eq. (2) where the dissipation kernel 7(i) 
depends on the bath oscillator frequencies and the cou- 
pling constants. Unlike the case with (5-correlated noise 
there exists no general proof that, for a general potential 
V{x), the system will reach thermal equilibrium at long 
times. The reason for this is that in this case the con- 
struction of a Fokker-Planck description is difficult and is 
known in few cases {e.g. harmonic oscillator case treated 
in Ref. [6]). For the special case of a harmonic poten- 
tial and with certain restrictions on the form of 7(i) one 
can prove equilibration by a direct solution of the equa- 
tions of motion [7] . In the quantum mechanical case the 
oscillator bath model has been widely used to model the 
effects of noise, dissipation and decoherence in quantum 
systems [8,9]. In this case the approach to equilibrium has 
been proved only for a special class of potentials for the 
cases where the system-reservoir coupling is weak [10]. 

A number of papers [11-13] have attempted to under- 
stand various aspects of the generalized Langevin equation 
such as anomalous diffusion, nonstationarity and ergodic- 
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ity. In this paper we address the question of approach to 
equihbriuni for the generahzed Langevin equation. For a 
particle in a harmonic potential V{x) — lUqx'^ /2 and cou- 
pled to a heat bath with a finite band-width (and hence 
long-time memory) we show that equilibration is not al- 
ways ensured and depends on the oscillator frequency luq. 
We find the necessary conditions for equilibration which is 
related to the existence of bound states (localized modes) 
of the coupled system-plus-bath. We note that Ref. [13] 
looks at the question of ergodicity which is similar to the 
question addressed here. One of their results is that the 
motion of a harmonic oscillator coupled to a general heat 
bath is ergodic. This is based on the fact that the coef- 
ficients in the Fokker-Planck equation asymptotically ap- 
proach constant values. Our work shows that this condi- 
tion does not ensure ergodicity. Next we consider a par- 
ticle in a nonlinear potential V{x) = ujqx'^/2 + ux'^/4 and 
coupled to a special bath, the so-called Rubin model [5,14]. 
The usual expectation would be that nonlinearity should 
help in equilibration. However surprisingly we find the 
contrary to be true. For values of t^o for which the lin- 
ear system (with u = 0) does equilibrate, we show that 
switching on the nonlinearity can lead to loss of equilibra- 
tion. We relate this to the formation of nonlinear localized 
modes which arc similar to the discrete breather modes 
found in nonlinear lattices [15,16]. 

Definition of model. — We consider the following 
Hamiltonian of the usual model of a system coupled to a 
bath of N oscillators: 



N 



H = H., 



p2 






Cq^X 

UJ 



(4) 



where {Xa, Pa} denotes degrees of freedom of the a**^ 
bath oscillator (a = 1,2,...N), lo^ is its frequency, and 
Cq is the strength of the couphng to the system. For 
dissipation it is necessary that the frequencies Wq have 
a continuous spectrum in the limit N -^ oo. We as- 
sume that at time t — to the system and bath are de- 
coupled ({cq} = ). The bath is in thermal equilibrium 
and described by the canonical distribution e^^^'' where 
Hb — YlioX Pal'^ + ^li^al'^ ] E^nd the system is in an arbi- 
trary initial state {x(io), p(io)}- The coupling is switched 
on at time t ~ t^. It can then be shown that eliminating 
the bath degrees of freedom leads to the following equation 
of motion for the system (for t > to): 



dV{x) 
dx 



■-f{0)x- 



;.-M^.,., 



-Vit), (5) 



with the dissipation kernel given by 7(f) = 
J2ai^a/^a)^os{uJat). In the limit to —> -00, using 
7(t -^ 00) = 0, we get eq. (2). The noise correlations 
satisfy the generalized FD relation in eq. (3) where the 
noise average (...) is obtained by averaging over the initial 
conditions of the bath. 



Equilibration in a harmonic potential. — We first 
consider the case with V(x) — luqx'^/2. For this case 
eq. (5) is a linear inhomogeneous equation whose solution 
is: 

x{t) = H{t-to)x{to) + G{t-to)p{to) 

+ I dt'G{t - t')7^{t') , (6) 

where H{t) and G{t) are the solutions of the homogeneous 
part of eq. (5) (with to — 0) for initial conditions -ff (0) = 
1, ij(0) == and G{Q) ^ 0, G(0) = 1 respectively For 
equilibration we require that, in the long time limit, the 
solution in eq. (6) should not depend on initial conditions. 
Thus a necessary condition is that G{t -^ 00) -^ 0. At 
large times it can be shown that H{t) -^ G{t) and hence 
H{t -^ 00) -^ also. Let us define the Green's function 
G^{t) = G{t)9{t). It is easy to see that G^ satisfies the 
equation of motion: 



G^ 



- w^'G^ 



+ / dt'-i+{t-t')G+{t')^5{t) 



where 7^(t) ~ j{t)9{t). Defining the Fourier transforms 
G+{uj) = JZodte"^'C!+{t) and j+{uj) = /~^ dte''^*7+(t) 
we find that G+(a;) is given by 



G+(^) 



-(lU^ — UJq) — iLU"f+(uj) 



(7) 



with 7^(w) — IbX^) + "^ui^) where the real and imag- 
inary parts arc given by: 'j'jfioj) — J2ai''^^a/'^^a)[^{^ + 
iVa) + S{lj - LJa)], -fiiuj) = Eaclu}/{^1(^^~^1)]- 
Now from eq. (7) it follows that G~^{t -^ 00) — 
limf^oo(l/27r) J(icL;e^*'^*G+(a;) will vanish (this follows 
from the Riemann-Lebesgue lemma [17]) unless G^{lo) 
blows up at some real value of w. Thus the condition 
under which G+(t -^ c>o) does not vanish is that the equa- 
tion 



-(uj — cJq) — iuj^ (bj) = 



(8) 



has a real uj solution, which we will denote by flf,. It 
follows from eq. (8) that "fji{Q,i,) ~ which means that Qb 
necessarily lies outside the bath band-width. This solution 
corresponds to a bound state. To see this we now solve 
the equations of motion of the coupled system-plus-bath 
by a resolution into its normal modes. Let us denote the 
normal mode frequencies by fig and the corresponding 
eigenfunction by Ug == {Uo^Q,Ua=i^Q, ■■■Ua=N,Q} where 
Uqq corresponds to the system variable x. They satisfy 
the equations: 



a ^ 

CaUoQ 



^n%a 



Q'^aQ 



-LOI [ U. 



aQ 



< ) 

a^ 1...N (9) 
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A bound state [18] occurs if there is a mode such that 
its frequency Qt lies outside the bandwidth of the iso- 
lated bath. The corresponding eigenfunction Ut, will have 
a finite weight at the system point {i.e. Uot = 0(1)). 
Solving, for Uab, the second equation in eq. (9) and sub- 
stituting into the first equation, we find that the condi- 
tion for a bound state is given precisely by the solution 
(if it exists) of eq. (8). The corresponding eigenfunction 
is given by: Uab = CaUah/{-i^l + ujD, a = l...iV and 
Uob = [1 + Eccl/i-^l+^DT^^^- The general solu- 
tion of the full equations of motion in terms of normal 
modes is 



Y ^ Ucosn{t - to)U-^Y{ta) + U 



-IW. ^ ^rr^^^('-to\u-'Y{to) 



n 



where Y = {x,Xi, ...Xa...XN}"'" . From this one can iden- 
tify the part of x{t) involving p{to). Comparing with 
eq. (6) we get G{t) ~ ^g UQQsin{nQt)/QQ. If there are 
no bound states then the frequencies Qq form a continu- 
ous spectrum, the sum can be converted into an integral, 
and in the limit t ^ oo the infinite oscillations lead to the 
integral vanishing. In the presence of a bound state we 
get a non- vanishing contribution given by 



G{t 



2 sin(Ohi) 



-'06" 



n. 



(10) 



Thus we again arrive at the conclusion that in the presence 
of bound states, thermal equilibration is not achieved and 
the steady state properties depend on the initial conditions 
of the system. The long time form in eq. (10) can also 
be obtained directly from eq. (7) by integrating around 
the singularity given by eq. (8). In the absence of bound 
states, eq. 6 gives the unique steady state solution x{t) = 
/_ dt'G^{t — t')ri{t'), and one can verify that: 



fceT 



{^)^^rd.\G-i.)\^Re[,-i.)] . ^ 



^2' 27r 



duj oj^\G+{uj)YReY^+{uj)] = 



knT 



where the integrals are easy to evaluate if one chooses ap- 
propriate contours [7]. The inset in fig. (1) gives a numer- 
ical demonstration of the nonequilibration problem (see 
following section for details of numerics) . 

Equilibration in a nonlinear potential. Usually 
one expects that introducing nonlincarity should help in 
equilibration and we will now test this expectation. We 
will consider a particular model of a heat bath for which 
it will be possible to analyze things in detail. This is 
the Rubin model in which the bath is a one-dimensional 
chain of coupled oscillators. Thus the full Hamiltonian 
is: H^H, + Eti b?/2 + (a;.+i - X,) 2/2] + {x^ - xf 12, 
with xn+1 = and where x ~ Xf) refers to the system. 
Transforming to normal mode coordinates of the bath we 
recover the form eq. (4) with uja = ^^{^^{qa/2) and Cq = 
[2/{N + l)]i/2sin(qa) with qa = aT:/{N + 1) and a = 



1, 2...N. From this we can find the explicit form of j^{lu) 
which is: 

iujj+{uj) = -{l-e:'>) for |w| = 2|sin(g/2)| < 2 
= -{1 + 6-") for |a)| = 2cosh(cr/2) > 2. (11) 

For a harmonic potential V{x) — ujqx'^/2 we see from 
cq. (8) that the condition for getting bound states is 
luq > 2. The bound-state mode is given by x^ — 

A(— l)*e~'^' cos(fi(,i)- where fit — tJo/V^o ^ 1 i^ the 
bound state frequency, cr = In {loq — 1) and A is an ar- 
bitrary constant which can be fixed by normalization. 

We now introduce the nonlinear term ux^/4 and study 
its effect numerically. Let us choose cjg = 1, in which case, 
in the absence of the nonlinear term, there are no bound 
states and we expect equilibration. In our simulation we 
consider a heat bath with A^ — 1000 particles whose ini- 
tial positions and momenta are chosen from the Boltz- 
mann distribution e^^''''^ with T = 1. The system has 
the initial state {p(0) — 0, x{0)}. We solve the equations of 
motion numerically with this initial state with the system- 
bath coupling turned on at time t = 0. We compute the 
expectation values Ke{t) = {p^) and Pe{t) = (x dH/dx) 
where (....) denotes an average over the bath initial con- 
ditions. In our simulations we averaged over 10^ initial 
conditions. We have checked that results do not change on 
increasing the bath size or number of aveages. For equili- 
bration we expect the steady state values K^. = P^ = T. In 
fig. 1 we see that for the initial condition a;(0) ~ 7.0 the lin- 
ear problem equilibrates as expected. On the other hand 
for a small value of the nonlincarity parameter u = 0.1 the 
system fails to equilibrate. For a different initial condition 
a;(0) =: 5.0 the system equilibrates. The loss of equilibra- 
tion infact occurs whenever x(0) is greater than a critical 
value Xc ~ 5.6. 

The dependence of equilibration on initial conditions 
can be traced to the formation of localized states that can 
be generated once we have nonlinearity. These nonlinear 
localized modes [15,16] are stable time-dependent solutions 
of the equations of motion that are localized and largely 
monochromatic and have the same form as the impurity 
bound states namely x\{t) = A(— l)*e^'^' cos(r2tt). If we 
plug this into the equations of motion and neglect terms 
containing higher harmonics (^ Sfif,) we then get the fol- 
lowing conditions relating fib, a and A (for the parameter 
choice ujf) = 1): 



A2 = ^, VLl = 2- 
■iu 



(12) 



Note that, of the three parameters needed to describe a 
bound state namely flb,cr,A only two get fixed and we 
thus have a continuum set of possible bound states which 
can be described by a single parameter. It is clear that 
for any non-zero u we can always choose cr > and deter- 
mine A, rif, from eq. (12). Hence bound states exist for 
any value of u. A clean way to numerically observe these 
localized modes is by initially preparing the bath at T = 
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Fig. 1: Plot of average kinetic energy Ke{t) as a function of 
time for different initial conditions. The purely linear case 
always equilibrates (for Uo ~ 1) while in the nonlinear case, 
equilibration depends on initial conditions. The behaviour of 
Peit) is similar. The inset shows the nonequilibration prob- 
lem in the linear case. The parameter value ujq = 3.0 gives a 
bound state with fit — 3/2^" and we plot Ke{t) at large times 
for 2;(0) = 5.0 and T — 1. The solid line gives the analytic 
prediction K^{t) = G{tf x{Of +T with G{t) given by eq. (10). 



and the system with initial conditions {p(0) = 0,a::(0)}. 
Then at very long times we find that this initial con- 
dition relaxes to a localized mode if a;(0) is sufficiently 
large. In fig. 2 we plot the coordinates {x{t),Xi{t)} for 
i — 1,2,3,4 as functions of time (for parameter values 
(jjQ — 1 and u = 0.1) . For the initial condition a;(0) — 7.0 
it is clear that the long time solution is quite accurately 
described by the localized mode solution. In the fits with 
the localized mode we fix A from the numerically obtained 
amplitude of x{t) and obtain a and O;, from eq. (12). 
For the initial condition x{G) = 5.0 we find that the ini- 
tial energy quickly dissipates into the bath and at long 
times is distributed uniformly. A more accurate form for 
the localized mode can be obtained by including higher 
harmonics. At the next order, the solution is given by 
x\ = ^(i)(-l)^e-'^'"^cosf7bi-H^(2)(_i)ig-rT(^)icos3f76i, 
where again there is only one unknown constant. In our 
example we find that A^'^^ « A^^^ (For all values of u we 
get y4(2)/^(i) < 1/21). 

Thus we have shown that the finite-temperature non- 
equilibration problem is related to the formation of non- 
linear localized modes which, unlike in the harmonic case, 
depends on initial conditions. 

Discussion. — In summary we have examined the 
conditions necessary for thermal equilibration of a par- 
ticle in a potential well and evolving through a general- 
ized Langevin equation. For the case of a harmonic po- 
tential we show that for finite band-width baths the sys- 
tem will not always equilibrate. Using the microscopic 
model of a bath as a collection of oscillators we show that 
non-equilibration arises because of the formation of bound 
states in the coupled system-plus-bath. Surprisingly we 
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Fig. 2: Plots of the positions of various particles x{t) and Xi{t) 
for i = 1,2,3,4 as functions of time plotted after the system 
has reached a steady state. Initial condition was a;(0) = 7.0 
and all other positions and momenta set to zero. The solid 
lines correspond to the analytical prediction for the breather 
mode with A = 6.44, a = 1.135 and fit = 2.331. 



find that making the system nonlinear does not restore 
equilibration. On the contrary nonlinearity assists in the 
formation of localized modes and hence causes loss of equi- 
libration for arbitrary initial conditions. These nonlinear 
localized modes are usually studied in the context of peri- 
odic nonlinear lattices, while here they arise in a situation 
where only one spring is nonlinear. In this paper we have 
considered a special one-dimensional bath. However lo- 
calized modes also occur in higher dimensional nonlinear 
lattices with finite band-widhts [16] and hence we expect 
that the problem of equilibration is quite general. 

In Ref. [19] electronic systems described by the tight 
binding noninteracting Hamiltonian were discussed and it 
was shown that the formation of bound states with finite 
band width reservoirs leads to a similar problem of equi- 
libration. One example which was studied was that of 
a single site attached to a one-dimensional free electron 
reservoir. An interesting question that the present study 
raises is whether the introduction of an onsite interaction, 
say of the Hubbard type, would lead to similar effects as 
the nonlinear term in the case of the oscillator. 
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